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The non-crossing approximation (NCA) is generalized for the Anderson model with finite 
Coulomb repulsion U. Resummation of infinite number of terms is performed so as to incorpo- 
rate all leading contributions in the limit of large degeneracy of the local states. With negligible 
weight of the doubly occupied local states in equilibrium, the extension is achieved through 
simple modifications of the lowest order formulae. The single-particle excitation spectrum is 
calculated with almost the same computational effort as in the original NCA. The present 
scheme reproduces the scaling behavior of the Kondo resonance in the density of states, and 
gives reasonable description of thermodynamics of the finite-f/ Anderson model over a wide 
range of temperature. 
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1. Introduction 

Heavy fermion systems, where 4/ electrons play an 
important role, show plenty of fascinating phenomena. 
Such systems are often modeled by the periodic Ander- 
son model, which exhibits the Kondo effect as well as 
the mixed-valence property. The single-impurity Ander- 
son model is utilized in the intermediate stage for solving 
the periodic model by means of the dynamical mean field 
theory (DMFT).^ Combined with band calculations, the 
DMFT can give a semi-quantitative description of den- 
sity of states in real compounds. For this purpose, it is 
required to solve the impurity model under realistic con- 
ditions for system parameters, which may lead to a heavy 
load in numerical calculation. Hence it is still important 
to develop a simple but reliable computational scheme to 
solve the impurity Anderson model, although low-energy 
properties of the model have now been understood well. 

Perturbation theory with respect to the hybridization 
can take into account strong correlations between local 
electrons. The lowest order self-consistent approxima- 
tion is called the non-crossing approximation (NCA), 
since the scheme incorporates all diagrams without cross- 
ing of conduction lines. In the case of large degeneracy n 
of the local states, the NCA derives proper dynamics as 
well as thermodynamics of the Anderson model with in- 
finite Coulomb repulsion U. The NCA accordingly plays 
a role complementary to the numerical renormalization 
group (NRG)i° and the quantum Monte Carlo (QMC)," 
since the latter numerical methods are only available in 
the small n case by restriction of computational time. 
Combination of the NCA and the DMFT has so far re- 
vealed various characteristics of the periodic Anderson 
model^^'^^ and the Hubbard model.^'' 

The NCA is originally derived for the infinite-?/ Ander- 
son model, where only 4/° and 4/^ states are contained. 
For practical use, however, additional informations on 
4/^ states are necessary as well. Extension of the NCA 
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to the finite-L/ case has been tried in several ways."'^^"^^ 
It has been found that increase of the Kondo tempera- 
ture, which should occur due to the involvement of 4/^ 
states, cannot be reproduced by lower order modification 
of the NCA to the finite-?/ case, even though n is large. 
Therefore, higher order terms are relevant in such mod- 
els. 

Sakai et al. have pointed out importance of infinite 
number of skeleton diagrams in large-n limit, and has 
called their approximation NCA/^v.^^ The theory seems 
to yield a proper energy scale of the Kondo effect. How- 
ever, the resummation procedure for the single-particle 
excitation spectrum does not give a conserving approxi- 
mation. Haule et al. have developed a conserving approx- 
imation by selecting a set of generating functionals.^^ It 
is called symmetrized finite-t/ NCA (SUNCA), since the 
approximation incorporates excitations of both particles 
and holes on equal footing, and can accordingly repro- 
duce the symmetric feature of the single-particle excita- 
tion spectrum. Although the theory provides a proper 
energy scale in n — 2, actual computation is not easy 
for large n even for a simple structure of the conduction 
band. 

In order to be applied to realistic impurity models with 
complex 4/-shell and band structures, or to the DMFT 
as a solver of the single-impurity problem, the practical 
calculational scheme should be simple enough. In this pa- 
per, we modify the NCA/^v scheme so as to recover the 
conserving property in the limit of large n, and present 
results for the single-particle spectra. The computations 
require almost the same effort as in the original NCA. 
It will be shown that the theory gives a proper energy 
scale and characteristics of the Kondo effect under the 
large-n condition. In addition we check the accuracy of 
the theory in thermodynamics. 

This paper is organized as follows. Section 2 introduces 
a model that we will consider in this paper, and gives a 
brief description of the resolvent method. Formulations 
for the single-particle excitation are given in §3. In §4, 
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we present numerical results, and discuss their accuracy 
in several aspects. In the final section, we summarize 
characteristics of the present theory. 

2. Model and Resolvent Method 

We consider the Anderson model with an orbital de- 
generacy in the local states, which are called 4/ states 
hereafter. Assuming strong Hund's coupling, we distin- 
guish the 4/ states by a single quantum number a with 
n-fold degeneracy, e.g., the azimuthal component of the 
total angular momentum j = 5/2. Conduction electrons 
hybridize with 4/ electrons having the same index a. 
Then, the Hamiltonian takes the form 



ha. 



u 



E flfo.fi' U' + E Vkflc^c + h.c, (1) 

a^a' ha 

where Cka and /„ is annihilation operators for the con- 
duction electron and the 4/ electron, respectively. 

We work under the large-J7 condition, and accordingly 
restrict 4/ states to 4/°, 4/^ and 4/^. In order to take 
account of the strong correlation between 4/ electrons, 
we introduce the resolvent R^{z) for each 4/ state as 
follows: 



Rj{z) — [z — — 



(2) 



where 7 indicates either Af^, 4/^ or 4/^ states, and here- 
after we use corresponding indices 0, 1 or 2, respectively. 
ej is the unperturbed energy, i.e., eo = 0, ei = e/ and 
€2 = 2e/ -|- [/. The hybridization is to be included in the 
self-energy part T,j(z) by the self-consistent perturbation 
theory.^"® 

Once the resolvent is obtained, physical quantities are 
computed with use of them.^"® For example, the partition 
function Zf of 4/ part is given by 

-^^Y.R,{z), (3) 



^^■=/^2^^ 



where the contour C encircles all the singularities of 
R-yiz) counter-clockwise. The resolvent is utilized to cal- 
culate dynamics as well. 

Throughout this paper, we give formulations on the as- 
sumption that there is no splitting in 4/^ or Af^ levels for 
simplicity. Alternative formulations with such splittings 
are obtained straightforwardly by distinguishing resol- 
vents for each split level. 

3. Finite-C/ Perturbation Theory 

The excitations of the conduction holes always accom- 
pany a summation of 4/^ or 4/^ states, which produces 
a factor n or n — 1, respectively. Hence from viewpoint of 
the \/n expansion, the hole excitations are not in higher 
order in perturbation series, and proper account should 
be taken of consecutive occurrences of such excitations. 
Since 4/°-4/^ fluctuation terms are represented by the 
diagrams without crossing of conduction lines, they are 
counted by the self-consistent treatment of the lowest 
order skeleton diagrams. On the other hand, 4/^-4/^ 
fluctuations involve crossing of conduction lines. Conse- 



z+e 



M e 
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z 



+ 



.-<--t^L^. 



Fig. 1. Diagrammatical representation of the equation for the 
vertex function. SoUd, wavy, dashed and double dashed hne rep- 
resent conduction electron, 4/*^, 4/^ and 4/^ resolvent, respec- 
tively. All resolvents are regarded as renormalized ones. 



quently, higher order skeleton diagrams are necessary to 
take account of such terms. 

In this section, we present an extended theory where 
4/^-4/^ fluctuations accompanying conduction holes are 
included properly. We first review the theory by Sakai et 
al. where a vertex part is introduced for the self-energy. 
After a brief description is given of a relation to the the- 
ory by Haule et at, we proceed to the formulation of the 
single-particle excitation. 

3.1 The NCAf^v scheme 

Before we present formulations for the single-particle 
excitation, we summarize the NCA/^v scheme developed 
by Sakai et al.'^' In order to derive the self-energy, they 
have introduced the vertex part A(z, e) which consists of 
the Af^-Af^ fluctuations with consecutive excitations of 
conduction holes: 

A(z, e) = 1 + (71 - 1) j dt'W{e')f{e')Ri{z + e') 

X R2{z + e + (')K{z,e'), (4) 

where /(e) is the Fermi distribution function and W{e) 
is defined by W{e) = J2k |Vfcp(5(e - efc). Equation (4) is 
diagrammatically represented by Fig. 1. The vertex part 
can equally be represented by diagrams with all arrows 
reversed in Fig. 1. 

A computation of the function A(z, e) with double ar- 
gument takes a lot of time. Sakai et al. have pointed 
out that the energy dependence of R2{z) is less impor- 
tant than i?i(z)'s in eq. (4), and proposed an approx- 
imation where R2{z) is replaced by a constant value 

i?2 = R2^\^f)i which is the bare resolvent of 4/^ state at 
z = e/ . In this approximation, the vertex part is indepen- 
dent of the energy of conduction electrons, and becomes a 
function of single argument written as A(z). Equation (4) 
is then solvable and leads to an expression 



A(z) = [l-(n-l)i?2S(^''^^(2 
where Sq^*^^'' (z) is defined by 



v(NCA), 
^0 



(z)= / deW{e)f{e)Ri{z + e), 



(5) 



(6) 



which is the self-energy of 4/" state in the NCA divided 
by n. The approximated vertex function A(z) is now an- 
alytic in the region except for Imz = 0, while original one 
A(z, e) has the additional cut at Im(z + e) = coming 
from i?2(z). 

The self-energy part of each state is represented with 
use of the vertex part introduced above. The self-energy 
of 4/° state is obtained by multiplying the NCA self- 
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(a) (b) 




Fig. 2. Diagrammatical representation of the self-energy (a) 
J2o{z), (b) T,i{z), (c) Tj2{z), and (d) a correction term which 
should be included in T:2{z) for construction of a conserving ap- 
proximation. 

energy and the vertex function A(z) as follows; 

^o{z)=nt'^^''^\z)A{z). (7) 

Corresponding diagrammatical representation is given by 
Fig. 2(a). Figure 2(b) shows a diagram of 4/^ self-energy, 
whose analytic expression is 

Si(^) = J deW{e)[l^ f{e)]Ro{z~e) 

+ (n -I) J deW{e)f{e)R2{z + e), (8) 

where Ro{z) is defined by 

Ro{z) = i?o(2)A(z)2. (9) 

Introduction of the auxiliary function Ro{z) makes 
eq. (8) look similar to the lowest order self-energy. The 
self-energy of 4/^ states are evaluated without any vertex 
corrections by 

£2(2) = 2 J deW{e)[l ~ f{e)]Ri{z - e), (10) 

which is expressed by the diagram of Fig. 2(c). The term 
shown in Fig. 2(d) is actually required for a construction 
of the conserving approximation as we will discuss later. 
We neglect the term, however, since its contribution is of 
order at most. Thus we avoid the double integral 

with respect to the energy of conduction electrons that 
would take much time. The self-energy in the NCA/^v 
consequently contains all diagrams of the leading con- 
tributions in the large-n limit, while the original NCA 
for the infinite- 1/ model incorporates diagrams up to the 
first order of l/n. 

Resultant NCA/^v equations have forms quite similar 
to that of the simplest finite-C/ extension in spite of the 
resummation of all the leading contributions. The only 
modification of the equations is the replacement of So(^) 
and Ro{z) by Eo(z)A(z) and Ro{z)A{z)'^ , respectively. 
The vertex function A(z) has been given analytically by 
eq. (5). Accordingly, the computational effort is almost 
the same as in the lowest order self-consistent approxi- 
mation without vertex corrections. 




Fig. 3. Eighth order diagrams of the generating functionals (a) 
considered in both the SUNCA and the NCA/^v, and (b) con- 
sidered in the SUNCA but not in the NCA/^v. Diagram (c) rep- 
resents one of contributions to the 4/^ state self-energy derived 
from the functional (b) , which is of order l/n^. 

3.2 Relation to the SUNCA 

Haulc and coworkers have extended the NCA to the 
finite-?/ case, and called the scheme SUNCA. They 
have constructed a conserving scheme by taking a set 
of infinite number of generating functionals into consid- 
eration. We discuss relation between the SUNCA and the 
NCA/^v. 

The generating functionals in the SUNCA consist of 
two parts except for terms less than the sixth order. Fig- 
ures 3(a) and (b) show diagrams of the eighth order gen- 
erating functionals that the SUNCA incorporates. The 
functional (a) generates the self-energy terms included in 
the NCA/^v. On the other hand, the NCA/^v does not 
incorporate contributions generated from the functional 
(b). Figure 3(c) shows one of the eighth order terms of 
the 4/^ state self-energy derived from the functional (b). 
It includes three conduction electrons propagating for- 
ward, and accordingly is of order Corresponding 
vertex part coming from the (b)-type generating func- 
tionals consists of the 4/°-4/^ fiuctuations with consec- 
utive excitations of forward propagating conduction elec- 
trons. Hence perturbation series of the self-energy gen- 
erated from the (b)-type functionals are of order l/n at 
most. Consequently, large-n condition justifies the ne- 
glect of (b)-type diagrams, which are relevant in the 
small-n case, especially n — 2. 

It is noticed that the NCA/^v itself is a conserving 
scheme provided no further approximation is made such 
as simplification of the vertex part and neglect of the 
correction term for T,2{z). Since these approximations 
are justified under the condition of large-n and large-C/, 
we expect that violation of the conserving property is 
not serious under such condition. 

3.3 Single-particle Green function 

We now proceed to the single-particle excitation. The 
4/ electron Green function is represented in terms of the 
X-operator X^j/ = |7)(7'| as follows: 



Gfir) = - E (7l/a|7')(7"l/Il7"') 

77'7"7"' 
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Fig. 4. Diagrammatical representation of the 4/ electron Green 
function. 

X (T,X^y(T)Xy.y"), (11) 

where X'y^i{T) is in the Heisenberg picture. The Green 
function is now independent of quantum number a. After 
the Fourier transformation, the Green function G/(ie„) 
is to be derived with use of the resolvents, where e„ — 
(2n + 1)7tT is the Matsubara frequency. 

In order to construct a consistent approximation, in- 
finite terms should be considered in the single-particle 
Green function as in the self-energy. The original 
NCA/^v scheme, however, does not take proper account 
of the vertex correction. An appropriate formula of the 
Green function is given by 

Gfiien) = ^[ ^e-f''Mz)R,{z + ie^) 
Jc ^"^1 

+ in-l)^ f ^e~P^R,{z)R,{z + ie^), (12) 
Zf Jc 27ri 

where the contour C encircles Imz = and Imz = — ie^ 
counter-clockwise, and Rq{z) is defined by eq. (9). Since 
the vertex function A(z) is independent of the energy of 
the conduction electron and has the same analyticity as 
Ro{z) in the present approximation, the Green function 
has been given simply by replacing Ro{z) with Rq{z) in 
the lowest order formula. Corresponding diagram is rep- 
resented by Fig. 4. We note that the first term of right- 
hand side of eq. (12) contains not only diagonal compo- 
nents {TrXQi{T)Xio) and {Tr X 12 {t)X 21), but also off- 
diagonal components {TrXi2{T)Xio) and {TrXoi(T)X2i) 
when the vertex part is expanded. 

Performing the integral with respect to z and ana- 
lytic continuation ie„ lu + iS onto real frequencies in 
eq. (12), we obtain the single-particle excitation spec- 
trum Pf{w) — —TT~^lmG f{uj + iS) as follows: 

Pf{uj) = J de[|o(e)?7i(e + t^) + ?7o(e)Ci(e + ^^)] 

+ {n-l) J de[ei(e)r?2(e + c^) + ?7i(e)6(e + ^)]- (13) 

We have included spectral functions of the resolvent 
77^ (w) and those of the defect propagator ^^(w) defined 

by 

r]^{uj) = -Tr^HmRj{uj + iS), (14) 

^-yiuj) = ZJ^e-'^'^T]y{Lu + iS). (15) 

Analogously, fjoiuj) and ^o('^) £^re defined from Ro{z). 
Since it is awkward to compute ^^(ct;) by eq. (15) due 
to the Boltzmann factor, £,-y{uj) should be calculated by 
another set of equations which are obtained by trans- 
forming the original ones (see Appendix). 
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3.4 Sum-rules 

The single-particle spectrum obtained within conserv- 
ing scheme satisfies the following exact sum-rules: 

J du;fiLo)pfiLo) - (Xn) + (n - 1)(X22> = ^, (16) 

J dLo[l ~ f{co)]pf{cu) = (Xoo) + {n - l)(Xn), (17) 

where rif is the average number of 4/ electrons. (X^^) is 
an average occupation of 7 state, and is given in terms 
of ^^(w) by 

(^77) = / ^-fi^)- (18) 

The sum-rule of eq. (16) follows from the fact that the 
left-hand side corresponds to G'/(— 0) in eq. (11). Simi- 
larly the left-hand side of eq. (17) corresponds to Gf{+0). 
Equation (16) has an intuitive interpretation that the 4/^ 
states are occupied in proportion to n{Xii), and the 4/^ 
states in proportion to n{n — l)(A'22)/2. Equations (16) 
and (17) lead to another sum- rule of the form 

J diupfiu;) = {Xoo) + n{Xu) + {n - 1){X22). (19) 

The sum-rule for £,-y{uj), given by eq. (A-8), ensures that 
the integral of Pfij-o) is always less than unity. It is clear 
that the lowest order formula of the Green function sat- 
isfies the sum-rule, eqs. (16), (17) and (19), as in the 
original NCA.^ The fact may imply that any correction 
terms in the Green function should not affect the inte- 
grated spectral intensity. 

We now discuss whether the density of states given by 
eq. (13) satisfies the exact sum-rule in the present theory. 
We can see from eqs. (5) and (9) that Ra{z) behaves at 
infinity as lim|^|^oo Ro{z) — z^^ , similarly to R^(z). The 
asymptotic behavior leads to the following sum-rule for 
Ro{z): 

J dLu ijoiio) ^ 1. (20) 

Equation (20) ensures that the particle number is not 
affected by the vertex correction, and is given by eq. (16). 
On the other hand, the approximate vertex correction 
changes the total intensity. Equation (17) is modified due 
to the vertex correction by 

J duj[l - f{Lo)]pf{Lo) = (loo) + (n - l)(Xn), (21) 

where (Xqq) is defined with use of £,Q{e) in a manner 
similar to eq. (18). Equation (21) follows from eq. (13). 
If (Xqo) were the same as (^oo), the exact sum-rule in 
eqs. (16), (17) and (19) would be satisfied. The spectral 
intensity in our framework is examined by numerical cal- 
culation in the next section. 

4. Numerical Results and Discussions 

In this section, we present numerical results for single- 
particle excitation spectrum and thermodynamic quanti- 
ties. Then we shall examine accuracy of the present the- 
ory. We take a rectangular model W{e) = Wq 6{D — \£\) 
with half width of conduction band being D. The degen- 
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Fig. 5. Single-particle spectra computed in the several schemes 
with parameters n = 6, Wo = 30, e/ = —1200 and U = 3|e/| = 
3600 at temperatures (a) T = 20 and (b) T = 0.1. 

eracy n takes n = 6 and n = 4. We take a unit such that 
D = throughout this paper. One may imagine the 
energy unit as of the order of Kelvin to meV for corre- 
spondence with actual systems. Parameters Wq, e/ and 
U are varied in the range where the 4/ particle number 
is almost unity. 

4-1 Single-particle excitation spectrum 

We first show numerical results for density of states 
Pf{io) of 4/ electrons. In addition to results in the present 
scheme, results in the other schemes are shown for com- 
parison: (i) p^p, which considers the vertex correction 
in the self-energy of resolvents but not in the Green 
function,^' (ii) NCA/^, which is the lowest order self- 
consistent approximation including 4/^ states but with- 
out the vertex correction, and (iii) infinite-t/ NCA, which 
is the original NCA with infinite Coulomb repulsion. 

Figure 5 shows single-particle excitation with param- 
eters n = 6, Wo = 30, ef = -1200 and U = 3|e/| = 3600 
at temperatures T = 20 and T = 0.1. In the NCA/^, the 
peak of 4/^ states is slightly sharper and is located at 
lower energy than the other schemes. There is no signif- 
icant difference between the present theory and p^ so 
far as large ui is concerned. Near the Fermi level, on the 
other hand, the heights of the Kondo resonance arc con- 
siderably different from each other as shown in the inset 



Fig. 6. Single-particle spectra near the Fermi level for several val- 
ues of U with n = 6, Wo = 30 and e/ = -1200 at T = 0.1. 

of Fig. 5(a). This is mainly due to different values for 

the Kondo temperature Tk in different approximation 
schemes. At much lower temperature T = 0.1, we can 

actually see from Fig. 5(b) that the heights of the Kondo 

(s) 

resonance peak are almost the same, except for Pj ■ It 
is known that the peak position with large n roughly 
gives Tk,^^ which also determines the zero temperature 
limit of the susceptibility and the Sommerfcld coefficient 
of the specific heat. The precise definition of Tk is dis- 
cussed in the next subsection. We notice that Tk in the 
present theory is larger than that in the infinite- 1/ NCA 
as we expect, while the energy scale becomes smaller in 
the NCAp. The inaccuracy of the NCA/^ scheme has 
been pointed out in ref. 4. 

In order to look through the Kondo behavior in this 
scheme, we show spectra near the Fermi level at low 
enough temperatures compared with the Kondo energy 
scale. Figure 6(a) shows spectra for several values of 
U with n = 6. In addition to results in the present 
scheme, p'^p is shown for comparison. Characteristic en- 
ergy scales turn out to increase as U decreases. We define 
the Kondo temperature Tk by the peak position of the 
resonance, and show rescaled spectra as a function of 
u)/Tk in Fig. 6(b). The spectra for U > 3|e/| are almost 
identical to each other, while the spectrum for f/ = 2|e/| 
is a little different from the others. The excellent scaling 
property of the Kondo resonance ensures an accuracy of 
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Fig. 7. Single-particle spectra near the Fermi level for several val- 
ues of V with n = 4, W^o = 45 and e/ = -1200 at T = 0.1. 



the calculated intensity in the present scheme. On the 
other hand, clearly fails to demonstrate the scaling 
property. The vertex correction in the Green function 
consequently turns out essential to reproduce the uni- 
versal behavior of the Kondo systems. Spectra of n = 4 
also exhibit the scaling property as shown in Fig. 7. 

Figure 8 shows Wq dependence of spectra with n = 
6 and U = 6000. The ordinate is the density of states 
multiplied by Wq. The peak hights in the scaled ordinate 
are nearly identical, which is also consistent with the 
scaling property as expected. 

4-2 The Kondo temperature 

We employ the definition of Tk as the energy scale 
of the system at zero temperature. Namely the 4/ 
magnetic susceptibility x at T = determines Tk by 



X = Cnf/TK, 



(22) 



where C is the Curie constant. Assuming n/ = 1 with 

frozen charge fluctuation, the Formi liquid theory gives 
the 4/ contribution to the specific heat coefficient 7 as 

^ = ^^7^ TT- 23 

61k n 

In the Fermi liquid theory, the quasi-particle density of 

states /9j(0) per spin-orbital channel at the Fermi level 
determines 7 as 7 = n7r^pJ(0)/3. On the other hand, the 
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Fig. 8. Single-particle spectra near the Fermi level for several val- 
ues of Wo at T = 0.1 with parameters n = 6, e/ = —1200 and 
U = 5|e/| = 6000. 



original density p/(0) of 4/ states is given exactly by the 
Friedel sum-rule as 

1 



P/(0) = 



7r2Wo 



^2 _ 



V n 



(24) 



The wave function renormalization factor a/ connects 
these densities of states by 



P/(0) = 



(25) 



At zero temperature, an approximate analytic formula 
for Pf{uj) has been proposed, which interpolates the 
Fermi liquid result and the high-energy tail of the Kondo 
resonance. Between the Fermi level and the peak position 
of the Kondo resonance, the density of states is approx- 
imately given by^®'^° 



P/M = — ■ 7 ,0 , Xo ' 



(26) 



where e/ and A are the renormalized 4/^ level and its 

width, respectively. Requiring the same number of 4/ 
electrons and 4/ quasi-particlcs, we obtain the relation 

AV(2/ + A2) = sin^ (7rn//n) . (27) 

Hence the renormalization factor is given by a/ = 
A/(7rWo). Furthermore comparison with eq. (23) relates 
the peak position if and Tk with n/ = 1 by^^ 

n^sin(7r/n) cos(7r/n) ^ 

Tr{n - 1) 

Hence we can derive Tk from the peak position of Pf{u>). 
The proportionality factor in eq. (28) is 0.99 for n = 6 
and 0.85 for n = 4. 

We next compare the Kondo temperature Tk derived 
from pf{u) with analytic formula derived by the scaling 
theory. In the parameter regime such that — gj 3> D 
and Cf + U ^ D, a combination of the Schrieffer- Wolff 
transformation and the second-order scaling theory gives 
an analytic expression for the Kondo temperature of the 
form^^ 
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Table I. The Kondo temperature for various sets of parameters. 
Tk is a value determined from the peak position of the Kondo 
resonance by eq. (28), and is that evaluated by eq. (29). 
(a) present theory (n = 6, Wo/\ef \ = 0.025) 













— GO 


2 


238 


292 


0.81 


— GO 


3 


77 


92 


0.84 


— GO 


4 


45 


52 


0.88 


— GO 


5 


33 


37 


0.90 


— GO 


oo 


10 


9.3 


1.09 


-1200 


2 


89 


- 


- 


-1200 


3 


23 


- 


- 


-1200 


4 


13 


- 


- 


-1200 


5 


10 


- 


- 


-1200 


oo 


5.5 


- 


- 


(b) prestnil i1kx..)i\' 


(11 = 


1 Mm / 


r r 1 O'-iT^^l 




U/\ef\ 


Tk 


Tk 


Tk/T^ 


— oo 


2 


201 


264 


0.76 


— oo 


3 


67 


80 


0.83 


— oo 


4 


40 


45 


0.89 


— oo 


5 


30 


32 


0.94 


— oo 


oo 


9.9 


7.9 


1.25 


-1200 


2 


70 




- 


-1200 


3 


18 




- 


-1200 


4 


10 






-1200 


5 


7.6 






-1200 


oo 


4.2 






(c) NCA/2 (n = 


= 6, Wo/\ef\ 


= 0.025) 




U/\ef\ 


Tk 


-'k 


Tk/T^^ 


— oo 


2 


7.9 


292 


0.027 


— oo 


3 


9.6 


92 


0.15 


— oo 


4 


9.8 


52 


0.19 


— oo 


5 


10 


37 


0.27 


— oo 


oo 


10 


9.3 


1.09 



g = nWo[\ef\-^ + {U + ef)-% (29) 

where c is a numerical constant of the order of unity, 
which cannot be determined by the scaling theory. We 
fix c as unity in this paper for simplicity. If the parame- 
ters are such that — e/ < D or tf + U < D, conduction 
band states with excitation energies around D do not 
participate in renormalization processes for the Kondo 
screening. Then the cut-off energies are given roughly by 
max{— for the lower and min{£),e/ + U} for the 
upper. In such a case, estimate of Tk becomes more com- 
plicated.^^' In order to avoid the complexity, we solve 
integral equations under the Coqblin-Schrieffer limit, 
that is, —ef,U,Wo oo with g fixed. Then eq. (29) is 
valid for comparison with the present theory. The limit 
is taken at the stage of integral equations by introduc- 
ing auxiliary functions such as Jo{z) = —WoRo{z) and 
110(2;) = So (.2)/ Wo that remain finite under the limit 
operation.^ 

(2) 

Table I shows Tk and ' , determined respectively 
by eqs. (28) and (29), for several values of parameters. 

Table 1(a) shows that the present theory gives a reason- 

(2) 

able agreement with in n = 6 over a wide range of U. 
Deviations from T^^^ become larger in the case of small 
U. We attribute the deviations to the approximation for 
the vertex function, which neglects energy dependence of 
the 4/^ resolvent in the renormalization processes. In the 



Table II. Integrated density of states of 4/ electrons in the 

present scheme evaluated by eqs. (16) and (21), and the exa<;t 
value by eq. (19), with ef = -1200 and T = 0.1. 



(a) n = 


6, Wo = 30 




U/\ef\ 


(present theory) 


(exact) 


2 


1.53 


0.962 


3 


1.13 


0.975 


4 


1.07 


0.983 


5 


1.05 


0.988 


(b)n = 


4, Wo = 45 




U/\ef\ 


(present theory) 


(exact) 


2 


1.48 


0.972 


3 


1.14 


0.983 


4 


1.08 


0.989 


5 


1.05 


0.992 



case of n = 4, on the other hand, variation of Tk/^k 
for different values of U/\ef\ is larger than the case of 
n = 6. This is mainly attributed to neglect of higher 
order terms in the 1/n expansion. Namely the present 
theory (U < 00) neglects many contributions of order 
1/n and beyond, while the original NCA ([/ = 00) keeps 
all 1 /n order terms but neglects some l/n^ contributions. 

We can see from Table 1(a) and (b) that the 4/*^ 
and 4/^ levels lying within the conduction band lead to 
quite lower values of Tk as compared to the case of the 
Coqblin-Schrieffer limit. The reduction of Tk should be 
due to the smaller cut-off energies caused by the virtual 
charge excitations. 

Table 1(c) shows Tk evaluated with use of the NCA/^. 
No enhancement of Tk by finite values of U is obtained by 
this lowest order self-consistent approximation. We thus 
conclude that inclusion of higher order terms of the ver- 
tex correction is the essence of the present theory which 
gives the proper energy scale of the Kondo effect. 

4-3 Sum-rules and exact relations 

We now examine whether the present theory satis- 
fies the exact sum-rules. Table II shows the integrated 
density of states in the present theory determined by 
eqs. (16) and (21), together with the exact value given 
by eq. (19). Note that the deviation from the exact value 
becomes larger as U decreases, while it is not affected by 
the value of n. We accordingly attribute the inaccuracy 
of integrated density of states to the approximation of 
the vertex function, because the NCA/^v scheme has a 
conserving property provided the integral equation (4) is 
solved exactly. We further conclude that the large inten- 
sity comes from the 4/^ part at high energy regime above 
the Fermi level. The conclusion is based on the following 
two reasons: first the particle number is not affected by 
the vertex correction, and secondly the Kondo resonance 
peak near the Fermi level shows the proper scaling be- 
havior. We have confirmed that the deviation from the 
sum rule becomes smaller at higher temperatures. 

Table III shows pf{0) obtained in the present scheme 
at T = 0.1 and the exact value at T = evaluated from 
nj by eq. (24) . The deviations from the exact values are 
roughly 10 percent in n = 6, and 15 percent in n = 4. 
The deviations are insensitive to the parameter U. Thus 
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Table III. The value pf{0)Wo obtained in the present theory at 
T = 0.1, and the exact value at T = evaluated from Uf by 
eq. (24), with parameter cj = —1200. 

(a) n = 6, Wo = 30 



U/\ef\ 


(present theory) 


(exact) 


Uf 


2 


0.0192 


0.0220 


0.926 


3 


0.0196 


0.0228 


0.944 


4 


0.0197 


0.0228 


0.944 


5 


0.0198 


0.0219 


0.942 


oo 


0.0198 


0.0219 


0.925 


(b)n = 


4, Wo = 45 






U/\ef\ 


(present theory) 


(exact) 


n f 


2 


0.0388 


0.0448 


0.926 


3 


0.0391 


0.0457 


0.938 


4 


0.0392 


0.0456 


0.936 


5 


0.0392 


0.0453 


0.933 


oo 


0.0391 


0.0441 


0.917 



the present theory is not worse than the original NCA 
in deriving p/(0). It is known that a fictitious anomaly 
appears at the Fermi level due to the inaccuracy of the 
NCA at r = 0.^"^ Wc can recognize tiny spikes at the 
Fermi level in Figs. 6 and 7. Note that the anomaly in 
n = 6 is smaller than that in n = 4. As in the original 
NCA, the present theory does not give exact description 
of the Fermi liquid property. 

4..4 Thermodynamics 

We now present numerical results for thermodynamic 
quantities of 4/ contribution, and then examine ade- 
quacy of the present theory for the Kondo effect. The 
partition function Zf of 4/ part is given in terms of the 
resolvents by eq. (3). The entropy Sf and specific heat 
Cf of 4/ contribution can be computed by numerical dif- 
ferentiations of Zf with resepct to temperature. To ob- 
tain reasonable numerical accuracy, however, it requires 
computations at numerous values of temperature, espe- 
cially for the second derivatives. Following ref. 24, we de- 
rive temperature derivatives of resolvents by solving ad- 
ditional integral equations. In this way thermodynamic 
quantities are computed without numerical differentia- 
tions. 

Figure 9(a) shows the entropy as a function of temper- 
ature for several values of U with n = 6. The entropy ex- 
hibits reduction from In 6 toward lower temperature due 
to the Kondo screening. The characteristic energy scale 
becomes large with decreasing U as in the single-particle 
excitation spectra. At low temperatures, the entropy is 
slightly less than zero because of the inaccuracy of the 
NCA and the present theory. On the other hand, values 
larger than In 6 at high temperatures come from charge 
fluctuations involving 4/" and 4/^ states. Figures 9(b) 
and (c) show temperature dependence of the specific heat 
Cf and Cf/T, respectively. The T-linear behavior of Cf, 
namely Cj/T being constant, is clearly seen at low tem- 
peratures. When temperature is scaled by Tk listed in 
Table I, C//T's with different values of U become almost 
identical to each other as shown in the inset of Fig. 9(c). 
Equation (23) gives ^Tk — 2.74 for n = 6, and is consis- 
tent with numerical results in the the inset of Fig 9(c). 
Thus the present theory with n = 6 demonstrates the 




0.1 1 10 100 1000 

T 




0.1 1 10 100 1000 

T 




' — 1 

0.1 1 10 100 1000 

T 

Fig. 9. Temperature dependence of the entropy Sf and the spe- 
cific heat Cf of 4/ contribution for several values of U with 
n = 6. Other pajrameters are same as Fig. 6. The values of Tk 
axe listed in Table 1(a). 

scaling property of the Kondo effect in the specific heat 
as well as the single-particle spectrimi. 

On the other hand, thermodynamic quantities with 
n = 4 show larger inaccuracy. The minus values of the 
entropy at low temperatures are seen more clearly than 
the n = 6 case in Fig. 10(a). Although the scaling behav- 
ior retains as shown in the inset of Fig. 10(c), Cf/T does 
not stay constant with decreasing temperature. Thus in 
the case of n = 4, the present theory (and the original 
NCA) has a limited accuracy to describe the specific heat 
at low temperatures. 
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Fig. 10. Temperature dependence of the entropy Sf and the spe- 
cific heat Cf of Af contribution for several values of U with 
n = 4. Other parameters are same as Fig. 7. The values of Tk 
are listed in Table 1(b). 



Wc obtain the relation between x and 7 by combin- 
ing eqs. (22) and (23) in the form independent of Tk as 
follows: 



X 



1 



(30) 



Derivation of the magnetic susceptibility for U < 00 \s 
much more complicated than the U = 00 case because of 
the 4/^ contribution. Hence we check the accuracy in de- 
scription of thermodynamics by returning to the original 
NCA {U = 00). Figure 11 shows temperature depen- 
dence of the susceptibility computed in the NCA. Prom 
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0.18 
0.16 
0.14 
0.12 

0.1 
0.08 
0.06 
0.04 
0.02 




n=6 



e^=-1200 
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10 
T 



100 



1000 



Fig. 11. Temperature dependence of the static magnetic suscep- 
tibility computed in the NCA with n = 6 and Wo/|e/| = 0.025. 
C is the Curie constant. 



Figs. 9 and 11 we obtain jC/x = 2.77 at T = 0.1, while 

eq. (30) gives 2.74. Thus the NCA {U = oc) describes 
with a reasonable accuracy the universal behavior in the 
specific heat and the static susceptibility. 

5. Summary 

We have derived the single-particle spectrum within 

the conserving scheme, following the NCA/^v approx- 
imation. The theory incorporates infinite skeleton di- 
agrams of leading contributions in the 1/n expansion. 
Equations have been simplified by approximations that 
are justifiable under the conditions of large-n and large- 
U. Resultant formulations have forms quite similar to 
those of the lowest order approximation: the equations 
for the Green function as well as the resolvents are ob- 
tained from the lowest order formulae by only replacing 
Eo(z) by So(z)A(z), and Ro{z) by Ra{z)A{z)'^. The ver- 
tex function A.{z) has been given analytically. 

Numerical calculations have shown that the scaling be- 
havior of the Kondo resonance is well reproduced by the 
present theory. Thermodynamical quantities computed 
in the scheme are properly scaled by the characteristic 
energy as well. Furthermore, the Kondo temperature Tk 
has turned out to be consistent with analytic estimate 
for n = 6 over a wide range of the Coulomb repulsion U. 
Although the present scheme breaks the sum-rule of the 
spectrum due to the simplification in actual calculations, 
the inaccuracy in intensity comes only from the 4/^ peak 
at high energy regime. Thus the drawback does not affect 
the density of states in the vicinity of the Fermi level. In 
addition, the inaccuracy remains small as long as U is so 
large that the relation |ej| < ef + U holds. This is the 
case in most of Ce compounds. 

The present theory has the advantage that it hardly 
requires extra time for computations as compared to 
the lowest order self-consistent approximation. Conse- 
quently, the theory is easily applicable to the DMFT, 
and to models with complex 4/ levels and realistic band 
structures. 
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Appendix: Equations for the Defect Propagator 

Single-particle excitation has been represented in 

terms of two types of spectral fimctions T]-y (lu) and £,-y{(-o) 
defined by eqs. (14) and (15). Although ^-y(w) is related 
to r]j{u)) analytically, it is difficult to compute ^-yiui) di- 
rectly due to the Boltzmann factor. Following ref. 23, we 
transform original equations for the resolvent into equiv- 
alent but more convenient forms, and calculate ^^(w) by 
numerical iterations. 

In order to make expressions simple, we in- 
troduce an operator V defined by VRj{u)) = 
—Zj^c~^'^Tr~^ImRj{u} + iS). With use of the operator 
V, we define cr.y(w) = PT,j{u)) and ^a('^) = 'PA(u;). Op- 
erating V on eqs. (2) and (5), we obtain 



\R-y{LJ + i^)|^cr-y(a;). 



(A-1) 



a(u') - {n-l)R2\Aiu; + i5)\^a^^''^\u;), (A.2) 



where al^'"^\aj) is defined by aQ^^'^^uj) = VT,, 
and is calculated by 



-.(NCA), 



^(NCA) 




{CO) 



~(NCA) 



(UJ) 



deW{e)f{-e)^i{LO + e), 



(A-3) 



which follows from an operation of V on eq. (6). Similarly, 
operations of V on eqs. (7), (8) and (10) yield equations 
for a^{co) as follows: 



ai{oj) = J deW{e)[f{e)io{uJ-e) 
+ (n-l)/(-e)6(a; + e)], 



(72(0;) = 2 / deW{e)fie)Uij-e), 



(A-4) 



(A-5) 

(A-6) 



where £,o{uj) is defined by ^o{uj) = PR{){co). Equation (9) 
leads to an explicit expression 

|o(w) -Co(t^)Re[A(o. + 1,5)2] 

+ 2Re[Ro{u> + i5)]Re[A{u + i6)]^A{u)). (A-7) 



These linear equations, (A-l)-(A-7), do not determine 
norms of ^-yiu}), a^{co). The norms can be determined by 
the following sum-rule: 



J da;^^^(a;) = l. 



(A-8) 



This identity follows from eq. (3). Normalizing by 

eq. (A-8) is equivalent to determination of the partition 
function, and thus thermodynamics is obtained by a com- 
bination of r]j{uj) and ^^{co)?^ 
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